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Abstract 

The problem considered is the development of the necessary plate 
stiffnesses for use in the general purpose program VICONOPT for buckling 
and vibration of composite plate assemblies. The required stiffnesses 
inc ude the effects of transverse shear deformation and are for sinusoidal 
response along the plate length as required in VICONOPT The method is 
based on the exact solution of the plate differential equations for a 
composite laminate having fully populated A, B, and D stiffness matrices 
which leads to an ordinary differential equation of tenth order. 


Introduction 

The VICONOPT program is described in reference 1 and is summarized in 
reference 2. The program performs buckling and vibration analysis of any 
prismatic assembly of composite plates and, optionally, the optimization of 
such plate assemblies. The assumptions made in the analysis are that the 
response in the longitudinal or x direction is sinusoidal and that individual 
plates havfe stiffness properties that result from balanced symmetric 
laminates though asymmetric laminates are treated approximately For 
orthotropic plate assemblies without shear loading, the sinusoidal response 
is exact for simply supported end conditions. For other cases, a series of 
sinusoidal modes is used with a Lagrangian multiplier technique to obtain 
results for quite general support conditions and loadings. The use of the 
term "exact" in the title refers to the fact that the analysis uses stiffness 
matrices that result from the exact solution of the uncoupled inplane and 
out-of-plane plate differential equations derived from classical plate 
theory. Because each of these equations are of fourth order, analytical 
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expressions can be written for these stiffnesses. The introduction of 
transverse shear deformation increases the order of the out-of-plane 
differential equation to six so that a numerical approach is required to 
obtain its solution and the subsequent plate stiffnesses. Once the 
numerical approach is chosen, it can be developed in a general way to 
include the coupled case which is of tenth order when transverse shear is 
included and of eighth order otherwise. The coupled case allows 
treatment of any laminate exactly with an arbitrary location of its 

reference surface. 

The standard numerical approach for solution of the coupled plate 
equations with or without transverse shear is to write the plate 
equilibrium equations in terms of displacements, assuming a sinusoidal 
response in one coordinate direction to obtain a set of ordinary differential 
equations with independent variable the coordinate in the other direction. 
The characteristic roots are obtained by setting to zero the determinant 
the coefficients of the displacement variables. This approach leads to 
increasingly long and complicated expressions that are impractica to 
implement for the most general cases. For plates with fully populated , 

B and D stiffness matrices, the number of terms is of the order ten 
thousand when transverse shear is included. The approach herein is to 
write the differential equations as a first order system with the 
displacements and the associated forces as unknowns as was done by 
Cohen in reference 3. The resulting matrix has typically two or three 
terms in each matrix element for the most general case. These equations 
can be arranged so that the characteristic root of the differential equation 
system appears only on the diagonal, thus a standard linear eigenva ue 
solver can be used. The eigenvectors then give the relationships between 
all the forces and displacements. 

Because exact expressions are used for stiffnesses, there is no need to 
subdivide a plate into segments as is done with conventional finite e ement 
methods to obtain accuracy. Nodes are required only where different 
plates are joined, usually at an angle with one another. The assumptions of 
plate theory are such that continuity of rotation at the junction of two or 
more plates requires that the shear angle in the plane of the edge formi £8 
the junction must be zero unless all plates are coplanar. ey s ep in 

analysis is to use the shear angle as the additional unknown over classical 
platJ theory, as was done by Cohen in reference 3. After calculation of the 
plate stiffness matrix, the rows and columns corresponding to the shear 
angle are simply removed so that the remaining unknowns are the same as 
classical plate theory. Thus the part of the computer code involving 
assembly of the global stiffness matrix from individual plate stiffnesses is 


2 



unchanged from that used for classical plate theory. Only additional 
routines to calculate the plate stiffnesses including transverse shear need 
to be developed. 


Theory 


Governing Plate F.mmtionc 

It is assumed that the plate has fully populated A, B, and D stiffness 
matrices with a reference surface arbitrarily located in the x-y plane as 
shown in figure 1. The plate is loaded by uniform inplane stress resultants 
N x , Ny, and N X y that act in a centroidal plane shown in figure 1(a). The 
centroidal plane is located at a distance z c from the reference surface and 
passes through the centroid determined from beam theory considering the 
plate to be wide beam whose length is in the x direction. Figure 1(b) 
shows the additional inplane forces n x , n y , n xy , moments m x , m y , m X y, 

and transverse shearing forces q x , q y that occur during buckling or the 
amplitudes of these forces and moments when vibration at a frequency of 
(o is considered. The deflections u, v, and w of the reference surface are in 
the x, y, and z directions respectively. For vibration or buckling, deflection 
in all three directions is resisted by Winkler elastic foundations of stiffness 
Kx. Ky, and K z . The equations of equilibrium are 

n X’X + n xy»y * N x uq, xx - K x u + 4 7t 2 co 2 (mQ u - mi\j/ x ) = 0 

n xy>x + n y>y - N x vq, xx - K y v + 4 k ~ co 2 (mg v - mj\(/y) = 0 

<lx>x + q y ,y - K z w - N x w, xx - Ny w,yy - 2 N X y w, xy +4^ 0 2 m Q w = 0 

m xy>x + nr'y.y -qy - N x z c vg, X x + 4 rc 2 co 2 (mj v - m 2 \|/y) = 0 (1) 

m X’X + m xy’y 4x ' N x z c u(), xx + 4 ji 2 co 2 (m] u - m 2 y x ) = 0 

where a comma indicates differentiation with respect to the variables that 
follow, uo and vq are inplane deflections at the centroidal plane, mj is the 
jth moment of mass about the reference surface and vj/ x and \}/y are 
rotations of the normals to the reference surface about the y and’ x faces, 
respectively. These equations are generalized from those of reference 4 to 
allow for the applied load N x to act at the centroidal surface which may be 
different from the reference surface. As in reference 4, it is assumed that 
the only inplane force to affect the inplane equilibrium equations is N x . 

The inplane displacements in the centroidal plane, uq and vq, are given by 
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UQ = u ' Z C Vx 


( 2 ) 


VQ = V - z c fy 


Stiffness Matrix 

The objective of the analysis is to derive a stiffness matrix that relates the 
force quantities along the edges y = ± \ to the displacements along the 

same edges. (Because of the sinusoidal variation in the x direction, the 
stiffness matrix will involve the amplitude of these quantities). 

The desired displacement variables are 


d = 


i u 
v 
w 

V y 
i Yx 


( 3 ) 


where y x = w, x - Vx has been introduced as a fundamental displacement 
variable rather than the rotation Vx in order to preserve continuity of 
rotations along nodal lines corresponding to plate junctions. 

The corresponding force variables are 

"i n x y 

n y 

Qy (4 

my 

i m x y 


The imaginary number i has been introduced so that the phase shift 
between the various quantities that occurs for orthotropic plates without 
shear loading results in real plate stiffnesses tor such plates. I he 
transverse shearing edge force in the z direction, denoted as Qy, must be 
normal to the reference surface of the undeflected plate. This is 
accomplished by replacing q y with the Kirchhoff value 


4 



( 5 ) 


q y - Q y - m xy , x + Ny w, y + N xy w, x 


Note that the m x> ,, x term which is found in the Ktrchhoff shear term of 

classical plate theory is also present in the transverse shear case when v x 

is used as a fundamental displacement variable. (This fact can be shown 
from the principle of virtual work). 

The problem is changed to an ordinary differential equation in y by 
assuming a sinusoidal variation in the x direction. If the variables hi 
equations ( 3 ) and ( 4 ) are considered to be functions of y, the 
displacements and forces throughout the plate are given’ by 


Z(x,y) = exp 


1 71 




z (y) 


( 6 ) 


where 


y. = [f _ 


and X ts the hall-wavelength of the response in the x direction. The next 

lLl P f 1S l h °. CXpreSS aH unknowns in ter ms of z. The stress resultant strain 
relationships given in terms of the A, B, and D stiffness matrices are 

part, ally inverted to give needed quantities in terms of the fundamental 
“ive 0 /. l6rmS denVab ' e fr ° m fUndamemal ^.ble, without aTy y 


xy 


m 


‘xy 


h 1 1 h i 2 h ] 3 h i 4 h | 5 h | 6 

- h i 2 h 2 2 ^23 h 2 4 h 2 5 h 2 6 

" h 1 3 h 2 3 h 3 3 h 3 4 *135 1135 

^1 4 ' h 2 4 - h 3 4 h44 1145 l'4 6 

-h | 5 h 2 5 h 3 5 -h 4 5 h 5 5 h 5 6 
' h 16 h 2 6 h 3 6 -h4 £ h 5 5 h 5 5 


"xy 

K x 
m ,, 


m 


xy 


' q x ' 

h 7 7 h 7 8 

~Yx" 

- 7y/ 

" h 7 8 hg 8. 

.q y . 


( 7 ) 
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where 

e X = U ’X 

K x = "Vx»x 
Vx = w >x - Yx 


The constants iqj 


h;; appearing in the first of equations (7) can be calculated 
from ^he'sta'ndard A B, and D stiffness matrices of laminate theory^ The 
shear stiffness terms in the last two equations are determined from 

reference 5. 


I, is also necessary to write equation (5) without the appearance of any y 
derivatives. The result is 


fQ v - m,y„ + Ny (Sly - Yx l>78) + N * y) 
qy = 1 -Ny hg 8 


( 8 ) 


The five strain displacement equations are 


£ X y — V, x 

e y = v> y 
7y = w, y ' Vy 

Ky = - V y>y 

K X y = Yx.y ■ Yy.x ' 2 Vy.x 


(9) 


Using equations (7) and (8), the strain-displacement and 
equations can be written in terms of the elements of /. as 


equilibrium 

follows 


/' = P L 


( 10 ) 


where a prime denotes differentiation with respect to y. The elements of /. 
are assumed to be given by 


zj = cj exp 


P y N 


(ii) 


where B is a characteristic root of the differential equation with as many 
values as the order of the differential equation system. Substituting 
equation (11) into equation (10) results in 
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(R -pi)c=o 


( 12 ) 


eTumion nn e T e h n,ity r ,riX and c is ,he formed from the ct of 

are the eigenvalul^of T r °°' S ° f ' he differential equation 

multiplying dividing appropriate cl b I s C by ^ The' T * 

deformation' bases' 1 

easy to implement in a computer nmor- * ^ ^ not com Phcated and are 

associated eigenvector corresnnnrl P t S l™' F ° r 6ach ei § enva| ue, the 
eigenvector corresponds to the c; of equation (\ 1 \ if r • 

matrix with columns as J c H uanon (12). If C is a 

be associated with displacemen^n^ loL^half ^de’ f "° ted as a ’ wil1 
associated with forces as follows h f ’ denoted as b, will be 


■KJ 


LUJ 

Denoting quantities evaluated at v = -h n k„ , 3) 

evaluated at y = b/2 by superscript 2 Pip Y UperSCnpt 1 and quantities 

displacements and forces^, the t^ edge VZ olam ° f 'k 

b or the plate may be written as 

N 


d = 


I 

k=l 


a jk rk exp 


'I Pk 


N 


'j -^ a ji 


«* = 7,a jk r k expf^ 


k = l 

N 


f J ~2w b jk r 


( 14 ) 


k= 1 



N 


jk r k ex P 


r j pk 


k=l 
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. t be determined from the edge values and N is the 

where r^ are consta matrix form, equations (14) are 

order of the differential equation. In matrix form, 

written as 



d 1 " 

d2j 


'f 1 ' 

Lf2j 

Eliminating 


f 1 

_f 2 J 

where 

K i 


= E r 


= F r 


(15) 


(16) 


Eliminating r from equations (15) and (16) gives 


= K 


d 1 

Ld2 


(17) 


(18) 


The rows and columns of K -- S P;;t v 8 o :: ng Y ^H:^al e tn^o:t ^n as a for 
stiffness matrix of the same size an yie i ds a lower bound 

the classical case. An alternate P the ^ame final unknowns by 

solution is to reduce the .stiffness ™ ^ program; (For the 

setting m X y to zer0 w . . c U hstructuring capability of the 

case of two plates that are “P lanar ' of freedom at each edge to join the 
program may be used with five d g ^ is rem oved). As for the 

,wo plates before the extra J e f" e mmetr ' c for orthotropic plates without shear 
classical case, K is real and y resulting stiffness matrix can be 

^^(fuST^rrth e original coding for the classical case. 

Obtaining accurate numerical result^ for K f™' 1 |T “! thods use d to overcome 
problem for certain ranges of pafametera 
these problems is discussed in ppen 

F j p ^ r v a liies for Clampfd 

j \/irnMOPT described in reference 1 is 

•*, M a 
"etuencV"- vS ion' and noTTbe confused with the eigenvalues o, the 
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R matrix) that converge to the desired result For each trial v,i u 

finding the buckHnt or Vihr ? COmple, , e d,scussi °" »f 'he algorithm for 
obtain the eigen JSL ^ " W ° U ' d be quite difficult '» 

a m r g h neral c r by edg T e he for the 

exceeded * 3 '^it^^e^ ^ ^ are 

doubled to re, nr 8 to the original width, the P ”re oi" ri^Tcl'T 
determined " 11171 ^ ^ e ' genvalues exceeded for clamped edges to be 

determined as follows. 0 '' Ap^ndVTshL^ , 

is proportional to the plate width b, and thus the eigenvalues of* R are ' X 
propomonal to b. Noting tha, an eigenvalue equal to it corresponds to 

until 'alf the 1^,™"°" W ' lb Slmply su PP orled edges, successively halving b 

S » - 1 ” ■ »••• *• «« « . 

Results 

n.'zzx *• 

various boundary conditions literature tor 

sit 

is! - 

iength direction d ^ 

n. r eco' ^ d % f :rmr,r;he'"eri^r:L';-„rr-^ 

t c ty solutlon and the third column is the ratio of frequency n) to co cJ , 
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- frequency obtatned r"r T'.OS whTch h 

representative *of*»°relati»ely « tSTS^TjL 

still small until n becomes arg . _ ^ which is well beyond proportions 

small for the enure range. transverse shear is appreciable but the 

usually encountered, the effem two per cent. These results are 

present results have errors y reference 6 that first order shear 

in agreement with the cone » » / response such as 

deformation theory ts adequate P were obtained b y subdividing 

vibration or buckling. P nls and linking them together by 

the circular arcs into small g deerees of freedom at an edge. 

substructuring techniques retaining all five degrees 


Table 1 Accuracy of shear deformation 


theory for cylinder vibration. 



The problem considered is a flM panel wUh bl ^*® ^ web) simp i y supported 
designed to carry an axta o o, OOO^Ib m. ^ ^ has stiffness 

Wtde column of^ 30 mch^ng = ^ = 93xl0 6, E, 3 = E,2- 

ST» Bfl'a,. with units psk andjotsson rrnio v - 

klfm the ‘lamintne "nd' Subscript 3 refers to normal to the 



laminate plane The density is .0571 lb/cu in. The sandwich core material 
is available in a range of densities which result in a varying shear stiffness. 
A typical variation of the shear stiffness of a foam core suitable for 
sandwich construction is shown in figure 2. The foam core was assumed to 
have a constant Young's modulus of 1000 psi at all densities with Poisson 
ratio of zero. Using these properties, optimized designs were determined 
tor solid composite and sandwich construction assuming both classical and 
transverse shear deformation theory. Typical cross-sections resulting from 
these analyses are shown in figure 3 where the relative size of one bay of 
the wide column panel is shown along with the actual panel dimensions. 
Only the design based on transverse shear theory is shown for the solid 
composite as it differs little from that based on classical theory. The 
composite laminates considered were [ ±45°/0°] S for each plate with a 

central foam core added for the sandwich configurations. The thickness of 
the various layers are given using the subscripts 45 or 0 for the composite 
material and a c to indicate the total thickness of the foam core for the 
sandwich configurations. For sandwich construction, the efficiency is such 
that unrealistically high strains result if the design is unconstrained so 
allowable strain limits were imposed of .006 along a fiber, .004 normal to a 
fiber and .005 for shear strain. These limits led the optimizer to reduce 
t e thickness of the +45° plies to zero so they were omitted in the 
sandwich designs. 

The mass of the various designs as a function of core density is shown in 
figure 4. Designs based on classical theory show a monotonic increase of 
mass with core density since the shear modulus of the core does not affect 
the result. Accounting for transverse shear results in a minimum mass at 
about ten lb/cu ft core density. However there is a range of core densities 
that can achieve nearly the minimum mass. The upper line showing the 
mass of the solid composite panel is about 1.6 times the minimum mass 
achievable with sandwich construction. For the higher core densities and 
the solid composite, there is very little difference in the mass achievable 
using classical or transverse shear theory. However even for these cases, a 
panel that just meets the design buckling load using classical theory would 
have a buckling load 5 to 10 per cent below the design value if analyzed 
using transverse shear theory. 

The sandwich cross-sections shown in figure 3 were obtained for a core 
density of 6.9 lb/cu ft so they are in the range where transverse shear 
effects are significant. In this case, using classical theory to design the 
panel (the middle configuration of figure 3) results in a configuration that 



would buckle at only 65% of the buckling load design requirement if 
analyzed using transverse shear theory. 

frrll' Hin ? Remarks 

A numerical approach to obtain exact plate stiffnesses including ihe effects 
of transverse shear has been presented. The generality of the meth 
allows inclusion of fully populated A, B, and D stiffness matrices. These 
stiffnesses have been incorporated in the V.CONOPT Program for ana ys,s 
and optimization of composite panels. Before the incorporation of this new 
capability, VICONOPT was limited to plates analyzed with classical theory 
having a zero B matrix so the addition of the new capability is an 
improvement for cases when classical theory is used as well as the 

transverse shear case. 

Two key items in the approach are: 

(1) expressing the necessary equations as a first order system ,nv °^"S 
the displacement and force quantities desired in the final stiffness relatio . 

(2) using the shear strain rather than a rotation for the extra degree of 
freedom introduced by transverse shear. 

The first item allowed the characteristic roots of the differential equation 
to be obtained from a relatively simple matrix. This is in contrast to 
previous conventional approaches for this type of problem which could 
fead to matrices whose elements consist of over ten thousand terms for the 
tenth order system associated with the coupled case including transver ^ 
shear The second item allows the plate stiffness matrix to be reduced to 
the same size as that for classical plate theory because continuity of 
rotation requires the shear strain to be zero a the junction of two plates 
that are not coplanar. A junction having all plates coplanar is treated by 
creating substructures with all degrees of freedom present and eliminating 
the extra degrees of freedom before assembly in the final stiffness matrix. 
Thus no additional coding is required to assemble the global stiffness 
matrix and to solve the resulting eigenvalue problem. 

Results of the analysis applied to the vibration of a composite cylmder 
show that the shear deformation theory used gives results with little 
compared to the elasticity solution. The analysis applied to a typica 
design problem showed that sandwich construction can give significant 
mass savings compared to solid composite construction. The analyses 
makes possible the determination of optimum core density tor sandwich 



construction and indicates the need to account for shear deformations 
especially when the lower mass core materials are used. 
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Appendix A 

Matrix for Charact eristic Roots 

The eigenvalues of R of equation (12) are the characteristic roots of the 
differential equations governing plate behavior. The nonzero elements of 
the 10 by 10 matrix R are 


$ $ 

r l 1 = *66 = -« h l 3 


r ?3 


a a a z 

- r 86 = - 1715 = ' b 1 ! 06 = 'V h34 
r ll = A 6 = -bh2 3 
r 110 = r 56 = -bh3 6 


r ?3 = -Al = - ^ T 25 = 
r f 9 = -Al = -bh2 5 


a a z 

i71 07 = -lTh24 


r 33 = r 88 = 


r 35 = - r l 0 8 = 


_ «Nxyh88 


bh 78 


T tl = -r?9 = ahi 5 

43 = A 9 = *^ T 45 = 0 9 = ‘V h4 5 

rj 9 = bh 55 


T \2 = r ?6 = -a 

r 16 = -^33 

r ?9 = 46 = -bh3 5 
# # 

r 2 1 = r 67 = -ah 1 2 

r f 7 = -bh2 2 

r 2 1 0 = r 57 = -bh2 6 
b 

r 34 = r98 = § 

a 

r 3 10 = " r 58 = -^38 


r 410 = r 59 = bh 56 


ah 88 

S 


1 3 


r51 = r 6 10 = - ah l 6 

a 

r54 = -r9i0 = a + s 

a 2 h88 

r 510 = " bh 66 ' bS 

r?3 = 4l =^ r 65 = - b ri °l = 
r?2 = B y 

0 a 2 Nxyh88 
r 83 = B z + a X 2 + bS 

aN xyh7 8 

*85 = ’ r l 0 3 = ' abX 2 ' s 
bN v h78 

r 95 = -n 04 = - s 


N xv h88 

r 53 = * r 8 10 = b ( h 46 + S } 

b 7 8 . 

f55 = r l 0 10 = -a( h 46 + § ’ 


:?1 = -B, 


a 2 h 


1 1 


a^h \ 4 


+ cxXj 

rf 4 = -rfa = bXi 


«N XV 

r 84 = r 93 = s 

a 2 9 bN y 

r 94 = “^44 + b 2 X2+ s 

bN xv h ^ 8 

rj 0 5 = ' b2x 2 + b b 77 " § 


where 


a = 


7i b 

X 


S = 1 - Nyh88 
X\ = ti 2 (4 mi co 2 + 


Nx z c 

X 2 


) 


X 2 => I 2 t 


2 

4 m2 w 2 NvZr a 2 


+ bx2>^ h 4 4 


where 



Nv 

= 7t 2 b(4 mo to 2 +^2) ‘ bX p 
can be x, y or z. 


Equations for the classical case can be obtained by setting the transverse 
shear strains, y x and y y , equal to zero and noting that Vx - W ’X- The 
partially inverted stress strain relations corresponding to equations (7) are 
slightly different in that m X y and K X y are interchanged. Only t e irst our 
equilibrium equations from equations (1) are used (the fifth is satisfied by 
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incorporation into the final form of the third). Following the same steps as 
e transverse shear case, the matrix whose eigenvalues are the 
characteristic roots is of order eight. Elements with a superscript # 

caL V, !f Sy is 8 's V ub n ,rimed he f ‘ ranSVerSe / hear *PP'y also to the classical 
elements no. given above^are"'' ***' The remaini "S 


r l 4 = -f85 = -2ah3 5 
r 34 = r87 = b 

2q2 
"b~ hl 6 

r - _ M 2 « 3 

74 r 83 ~ a N X y - ^2 b46 


r 54 = " r 8 1 = - 


r 24 - -r86 = -2ah2 6 
r 44 = r 88 = 2ah5 6 

r 73 = B z + a 2 X 2 

o h 44-4h6 6 

r 84 = a 2 + b2x 2 + bN y 


Appendix R 

Numerical — Solution of Equation^ 

There are several possibilities of numerical problems arising in 
determining K given in equation (18) as 

K = F F‘ 1 

(Bl) 

First the eigenvectors of the matrix R which form the elements of E and F 

^“scaling Wdh ned aCCUrate , ly f0r f extreme ^portions and a wide range 
of scaling. With accurate values of the elements of E and F the solution 

procedure is as follows. Solve the system of equations 


E* x = F* 


The solution is 


x = (E*)-l F* = K* 


(B2) 


(B3) 


where superscript * indicates Hermitian transpose of the matrix The 
olution is accomplished by Gaussian elimination with pivoting on the 
argest current diagonal except as will be discussed subsequently for 
transverse shear cases. A typical element in row k in the matrix equation 

(B3) has a multiplying factor of exp( + =^) where Pk is the k«h characteristic 

root. If the imaginary part of Pk is too large, numeric overflow will occur 
To prevent this, a real number is subtracted from the argument of each 



exponential in the k>h row such that overflow will n01 occ “ r kth 

particular computer being used. This has the effect of multiplying 
equation by a constant which does not change the solution. 

An additional problem occurs for transverse becomes 

veTyTarge tnT the Te'Cderapproach the “jj™ ^J^nd 

Finite element formulations have had serious pro ®J" S f elemen(s 

having been accurately calculated for aluminum plates with width- 
thickness ratios as high as one thousand. 

In some cases, particularly for unloaded isotropic plates Cha ”“ 

roots P from equation 12 occur in repeated pairs and the above jrocedu 

ss rs rs “■= 

t rs g : ssrc 

from the perturbed foundation stiffnesses. 
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Figure 1. Positive direction of forces and moments per unit width act.ng on 
a plate element. 



Figure 2. 


Shear modulus variation of foam core material with density 


1 8 





Solid 




T 

1 


b 2 

_i 


I"* - b j ►! 

Sandwich 


Classical Theory 


] 


Sandwich 

Transverse Shear Theory 


Designs based on Classical Theory 


Configuration 



skin 


blade 


mass 


bj 

b 2 l 45 

l 0 

l c 

l 45 

l 0 

l c 

lb/sq ft 

Solid 

1.12 

1.43 .00518 

.0183 

- 

.00891 

.0526 

- 

1.199 

Sandwich 

2.43 

2.69 

.0091 

.1710 

- 

.0236 

.4488 

.606 



Designs Based 

on Transverse 

Shear Theory 



Configuration 



skin 


blade 


mass 


b 1 

b 2 l 45 

l 0 

l c 

l 45 

l 0 

l c 

lb/sq ft 

Solid 

1.11 

1.41 .00520 

.0179 

- 

.00868 

.0533 

- 

1.205 

Sandwich 

3.90 

3.07 

.0131 

.3914 

- 

.0226 

.6606 

.737 


Figure 3. Geometry of solid composite and sandwich configurations 
designed to carry axial loading of 5000 lb/in. Sandwich core density is 6.9 
lb/cu ft corresponding to a shear modulus of 7100 psi. Panel length is 30 
inches. All dimensions are in inches 
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